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Abstract 

We have searched for the decay B + — > ul + v in 78 fb _1 of T(4S I ) data (85.0 million BB events) 
accumulated with the Belle detector. The final state is fully reconstructed using the to decay into 
7T + 7T~7r and detector hermeticity to estimate the neutrino momentum. 155 ± 47 signal events are 
found in the data, corresponding to a branching fraction of (1.3 ± 0.4 ± 0.2 ± 0.3) • 10~ 4 , where the 
first two errors are statistical and systematic. The third error is due to the estimated form-factor 
uncertainty. (This result is preliminary.) 

PACS numbers: 13.25.Hw 
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INTRODUCTION 



The magnitude of V u b plays an important role in probing the unitarity of the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix At present, experimental constraints on \V u b\ come 
mainly from the analysis of the decay B — > X u lv. Both inclusive approaches (sensitive to 
all X u lv final states within a given region of phase space) and the exclusive reconstruction 
of specific final states have been attempted. For the latter. Belle has previously obtained 
preliminary results for the decays B — > ixlv and B — > plv |2[. This article presents a study 
of the decay B + — > ul + u, which has not been observed so far 0]. 

Using a 78 fb _1 dataset recorded on the T(4S) resonance (85.0 million BB events), events 
with a single (undetected) neutrino are selected and the neutrino momentum is inferred from 
detector hermeticity. The neutrino candidate is combined with an identified lepton (electron 
or muon) and an uj (reconstructed through uj — > 7r + 7r~7r°) and the final state is reconstructed. 
The B — > ujIv yield and the remaining background are determined by a binned maximum 
likelihood fit. 

EXPERIMENTAL PROCEDURE 
KEKB and the Belle detector 

Belle is located at KEKB, an asymmetric e + e~ collider operating at the center-of-mass 
energy of the T(4S") resonance. The detector is described in detail elsewhere It is a 
large-solid-angle magnetic spectrometer that consists of a three-layer silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cerenkov 
counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a 
super-conducting solenoid coil that provides a 1.5 T magnetic field. 

The responses of the ECL, CDC (dE/dx) and ACC detectors are combined to pro- 
vide clean electron identification. Muons are identified in the instrumented iron flux-return 
(KLM) located outside of the coil. Charged hadron identification relies on the information 
from the CDC, ACC and TOF sub-detectors. 

Dataset 

The T(4f>) dataset used for this study corresponds to an integrated luminosity of 
78.13 fb _1 and contains (85.0 ± 0.5) ■ 10 6 BB events. Continuum background is subtracted 
using 8.83 fb~ x of data taken below the resonance. 

A full detector simulation based on GEANT is applied to Monte Carlo events. This 
analysis uses background Monte Carlo samples, equivalent to about three times the inte- 
grated luminosity. Monte Carlo data for the signal decay, B + — ► ujl + v, are generated with 
three different form-factor models: ISGW2 (quark model UKQCD (quenched lattice 
QCD calculation 0]) and LCSR (light cone sum rules Q). To model the cross-feed from 
other B — ► X u lv decays, Monte Carlo samples generated with the ISGW2 and the De 
Fazio-Neubert model |8| are used. 
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Neutrino reconstruction 

Events passing the hadronic selection are required to contain a single lepton (electron 
or muon) with a cm. momentum p* greater than 1.3 GeV/c. In this momentum range, 
electrons (muons) are selected with an efficiency of 92% (89%) and a pion fake rate of 0.25% 
(1.4%). 

The missing four-momentum is computed for selected events, 

Pmiss = PHER + PLER Pi , 

i 

Emiss = E HER + E LER — ^Ei, (1) 

i 

where the sum runs over all reconstructed charged tracks and photons. The indices HER and 
LER refer to the high energy and the low energy rings, respectively. To reject events in which 
the missing momentum misrepresents the neutrino momentum, the following selections are 
applied. Events with a large charge imbalance are eliminated, \Qtot\ < 3e, and the direction 
of the missing momentum is required to lie within the ECL acceptance, 17° < 9 miss < 150°. 
The missing mass squared, rn 2 miss = E^ iss — p\ iss , is required to be consistent with the 
neutrino hypothesis, |m^ is J < 3 GeV 2 /c 4 . 

After applying these cuts, the resolution in p m i SS is around 140 MeV/c for generic B — > 
X u lv events. As the energy resolution is worse than the momentum resolution, the neutrino 
four-momentum is taken to be (p m iss, Pmiss)- The efficiency of the combined event selection 
and neutrino reconstruction cuts is about 17% for B — > XJ,v events. 



Final state reconstruction 

Pairs of 7's are combined to form 7r° candidates [E 1 > 30 MeV, 120 < 771(77) < 
150 MeV/c 2 ). The decay uj — > 7r + 7r~7r (branching ratio: (89.1 ± 0.7)% jpj) is reconstructed 
from all possible combinations of one 7r° with two oppositely charged tracks. Combinations 
with a charged track identified as a kaon are rejected and the following selections are im- 
posed: > 300 MeV/c, 703 < m(n + n~n°) < 863 MeV/c 2 . Combinations located far from 
the center of the Dalitz plot are removed by requiring the Dalitz amplitude, A oc \p n + x p n - | , 
to be larger than half of its maximum value. 

The lepton in the event is combined with the uj candidate and the neutrino, and the 
lepton momentum requirement is tightened, 1.8 < pf < 2.7 GeV/c. To reject combinations 
inconsistent with signal decay kinematics, the cut | cos 8 by I < 1-1 is imposed, 

2E* R Et -m 2 R -ml 

cos 9 BY = B Y - f Y - , 2 

2p B p Y 



where E B , p% and m# are fixed to their nominal values, VEherEler, y E b 2 — m B an d 
5.279 GeV/c 2 , respectively. The variables E Y , Py and my are the measured cm. energy, 
momentum and mass of the Y = uj + I system, respectively. For well-reconstructed signal 
events, cos Qby is the cosine of the angle between the B and the Y system and lies between 
— 1 and +1 while for background, a significant fraction is outside this interval. 
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For each ujIv candidate, the beam-constrained mass m& c and AE are calculated (E£ eam = 
V EherEler) , 

AE = E* beam - (E* + El + El) , (3) 

and candidates in the range > 5.23 GeV/c 2 and \AE\ < 1.08 GeV are selected. On 
average, 2.5 combinations per event satisfy all cuts and we choose the one with the largest 
LO momentum in the cm. frame. This choice is correct in 77% of the cases. 

Continuum suppression 

The background from continuum e + e~ — > qq events, q — u,d, s, c, is suppressed using 
three variables that exploit the fact that, in the T(45) frame, the two B mesons are produced 
nearly at rest and that therefore BB events are nearly spherical while continuum events have 
a more jet-like topology. These variables are: 

• The ratio R2 of the second to the zeroth Fox- Wolfram moment This ratio tends 
to be close to zero (unity) for spherical (jet-like) events. (The cut R 2 < 0.4 is imposed 
at the event selection level.) 

• The cosine of Othmst, where 9 t hrust is the angle between the thrust axis of the ul system 
and the thrust axis of the rest of the event. 

• A Fisher discriminant that selects events with an uniform energy distribution around 
the lepton direction jlH |. The input variables are the charged and neutral energy in 
nine cones of equal solid angle around the lepton momentum axis. 

A cut on the likelihood ratio combining the three variables is imposed. This selection is 
56% efficient for B + — > uol + v events while it eliminates 92% of the continuum background 
remaining after the R2 < 0.4 cut. 

The fit 

The signal yield and the remaining background are determined by a binned maximum 
likelihood fit in the AE vs. m(7r + 7r~7r°) plane taking into account the finite Monte Carlo 
statistics 9 bins in AE (bin width: 240 MeV) and 8 bins in m(7r + 7r~7r°) (bin width: 
20 MeV/ c 2 ) are used. This fit is performed simultaneously in three bins of lepton momentum, 
1.8 < pi < 2.1 GeV/c, 2.1 < y\ < 2.4 GeV/c and 2.4 < V \ < 2.7 GeV/c. 

Five components are fitted to the data: the B + — > ul + u signal, B —>■ X u lv background, 
B — ► X c lv background, fake and non-5 decay lepton background and continuum. The 
shapes of the first four components are determined from simulation, while the shape of 
the continuum component is given by the off-resonance data. The normalizations of the 
B + — > u)l + v, the B — * X u lv and the B — > XJu components are floated in the fit, while the 
other components are fixed (Table HI Fig. [TJ). 
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FIG. 1: The result of the fit assuming ISGW2 form-factors for B + —* uol + u. The AE 1 and 
m(7T + 7T~7r ) distributions in each bin are shown. The data points are continuum subtracted 
on-resonance data, the histograms are the components of the fit, as described in the text. 
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i.O < P; < z.i Liev /c 


z.i < Pi < z.4 (je\ / c 


Z.4 < < Z. f Lre V /C 


data 


1990 


667 


75 


B + -> ujI+v 


41 ± 13 


68 ± 21 


35 ±11 


B -> X u lv 


61 ± 28 


82 ±28 


21 ±5 


B -> X c Zi/ 


1743 ± 36 


415 ± 14 





fake, non .B 


19 ±3 


33 ±4 


3± 1 


continuum 


17 ± 12 


61 ±23 


9±9 


sum 


1881 ± 49 


659 ± 44 


68 ± 15 



TABLE I: The result of the fit assuming ISGW2 form- factors for B + — > ujI + u. For each p*[ bin, the 
number of events in the signal window, 763 < m(ir + ir^7r ) < 803 MeV/c 2 and \AE\ < 360 MeV, 
are shown for the data and the different components of the fit. 



form- factor model 


signal yield 


B{B+ -> ul+v) 


X 2 /ndf 


ISGW2 


144 ± 44 


(1.00 ±0.31) • 10~ 4 


1.05 


UKQCD 


145 ± 44 


(1.20 ±0.37) • 10" 4 


1.08 


LCSR 


176 ± 52 


(1.67 ±0.50) • 10" 4 


1.04 


average 


155 ± 47 ± 15 


(1.29 ±0.39 ±0.28) • 10" 4 





TABLE II: The yield in the signal window, 763 < m(vr + 7r^7r ) < 803 MeV/c 2 and \AE\ < 360 MeV, 
the corresponding branching fraction and the goodness of fit (estimated by the x 2 divided by the 
number of degrees of freedom). For each fit (using a given form-factor model), the error quoted 
on the signal yield and the branching fraction is statistical only. For the average, the first error is 
statistical and the second is the spread around the central value. 

RESULT AND SYSTEMATIC UNCERTAINTY 

The fit is repeated for each of the three form-factor models available for B + — ► ujI + v. 
For each model, the signal yield, N(B + — > ujI + u), is determined and the branching ratio, 
B(B + — > u)l + v), is calculated accordingly using the relation (Table lH|) 

N(B + -> wl + v) = N(B + ) x B(B + -> ul + v) x B(u -> tt+tt-tt ) x (e e + e M ) . (4) 

N(B + ) is the total number of charged B mesons in the data, B(u — > 7r + 7r~7r°) = (89.1 ± 
0.7)% and e e (e M ) is the mo del- dependent selection efficiency for ujtv (ujfiu) candidates. 
Averaging over the three models (giving equal weight to each), a branching fraction of 
(1.29 ± 0.39) • 10~ 4 is obtained. The spread around this average value, which amounts to 
0.28 • 10~ 4 , is used as an estimate of the form-factor model uncertainty. 

The quadratic sum of the experimental systematics (listed in Table IIII|) is 0.21 ■ 10~ 4 
or 16.4% of the branching fraction. The largest contribution is the uncertainty in the 
X u lv cross-feed. It is estimated by separately varying the fraction of B — > irlv and 
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value 


AB/B 


Ref. 


B° -> n'l+u 


(1.8 ±0.6) • 10" 4 


2.2% 


M 


5° p l + v 


(2.6lo; 7 ) • 10 


12.7% 


[9] 


other B — » X^z/ 




1.4% 




cross-feed (sum) 




13.0% 




neutrino reconstruction 




4% 




charged track finding (l,ir + ,ir~) 




3% 




cluster finding (tt°) 




4% 




(sum) 




9% 




XJv cross-feed 




2.8% 




lepton identification 




3.0% 




number of BB 


(85.0 ±0.5) • 10 6 


0.6% 




B{uJ -» TT+TT-TT ) 


(89.1 ±0.7)% 


0.8% 


[9] 


total systematic uncertainty 




16.4% 





TABLE III: Contributions to the systematic uncertainty. The size of each contribution is given as 
percentage of the branching ratio. The different components are discussed in the text. 

B — > pip decays (that are expected to dominate in the high pi region) within their re- 
spective experimental uncertainties. For the cross-feed from other B — ► X u lv decays, the 
fit is repeated modeling this component once with the ISGW2 (fully resonant B — > X u lv) 
and once with the De Fazio-Neubert model (fully non- resonant B — ► X u lv). Half of the 
difference between these two extreme cases is assigned as a systematic uncertainty. 

The next largest component is the uncertainty in the neutrino reconstruction, track find- 
ing and cluster finding efficiency. Other contributions taken into account in the calculation 
of the systematic error are: X c lv cross-feed (estimated by varying the fraction of B — ► D*lv 
in the X c lv component), lepton identification, the number of BB events and the uncertainty 
in the us — > tt+tt^tt branching fraction. 

CONCLUSION 

We have studied the decay B + — ► ul + v using 78 fb _1 of T(45) data (85.0 million BB 
events). The final state was reconstructed using the uj decay into it + -n~-n Q and detector 
hermeticity to infer the neutrino momentum. The signal yield and the remaining background 
were estimated by a binned maximum-likelihood fit. Repeating the fit for three different 
B + — > u)l + v form-factor models and averaging the result, 155 ± 47 signal events are found, 
corresponding to a branching fraction of (1.3 ± 0.4 ± 0.2 ± 0.3) • 10~ 4 , where the errors are 
statistical, systematic and estimated form-factor uncertainty, respectively. (This result is 
preliminary. ) 
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